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Abstract 

We give a extensive account of a recent new way of applying the Dirichlet form theory 
to random Poisson measures. The main application is to obtain existence of density for the 
laws of random functionals of Levy processes or solutions of stochastic differential equations 
with jumps. As in the Wiener case the Dirichlet form approach weakens significantly the 
regularity assumptions. The main novelty is an explicit formula for the gradient or for the 
"carre du champ" on the Poisson space called the lent particle formula because based on 
adding a new particle to the system, computing the derivative of the functional with respect 
to this new argument and taking back this particle before applying the Poisson measure. 

The article is expository in its first part and based on Bouleau-Denis |12| with several 
new examples, applications to multiple Poisson integrals are gathered in the last part which 
concerns the relation with the Fock space and some aspects of the second quantization. 

Keywords : Dirichlet form, Poisson random measure, Malliavin calculus, stochastic differen- 
tial equation, Poisson functional, energy image density, Levy processes, Levy measure, gradient, 
carre du champ. 

1 Introduction and framework. 

This lecture is an introduction to Dirichlet forms methods for studying regularity of random 
variables yielded by Levy processes, solutions of stochastic differential equations driven by Pois- 
son measures and multiple Poisson integrals. The main part of this study has been done in 
collaboration with Laurent Denis. 

A Dirichlet forms is a generalisation of the classical quadratic operator \Vf(x)\ 2 dx early 
introduced in potential theory. The concept has been developped especially by Beurling and Deny 
in the 1950's as an application of Hilbert space methods in potential theory, and by Fukushima 
in the 1970's in connection with symmetric Markov processes theory. It received recently a 
strong development in infinite dimensional spaces where it appears as an alternative approach 
to Malliavin calculus. 

The importance of the notion comes from the fact that if Pt is a symmetric strongly continuous 
contraction semigroup on a space L 2 (fi) (for fi a- finite positive measure) with generator A, 
a necessary and sufficient condition that Pt be Markov is that "contractions operate" on the 
quadratic form £[f] = — < Af,f >l 2 0) i- e - £[<^(/)] ^ £[f] f° r contraction from M to E (cf 
Chap. I prop. 3.2.1). Such a quadratic form is called a Dirichlet form. 



The case of Malliavin calculus is that of Wiener space taking for Pt the Ornstein-Uhlenbeck 
semi-group. The corresponding Dirichlet form 8 possesses a "carre du champ" operator, i.e. may 
be written £[f] = ^ fF[f]d[i where T is a quadratic operator from the domain of £ to 
This fact makes it possible the definition of a "gradient" satisfying the chain rule and allowing 
a differential calculus through stochastic expressions and stochastic differential equations (SDE) 
and providing integration by parts formulae which yield existence of density results (cf |29|). 

Using Dirichlet forms in this framework of Wiener space improves several results : contraction 
arguments show that the Picard iteration method for solving SDE's holds not only in L 2 but 
still for the stronger Dirichlet norm. This gives existence of density for solutions of SDE's under 
only Lipschitz assumptions on the coefficients (cf |10] and [H]). More generally, Dirichlet forms 
are easy to construct in the infinite dimensional frameworks encountered in probability theory 
and this yields a theory of errors propagation through the stochastic calculus (cf Bouleau [7]), 
also for numerical analysis of PDE and SPDE (cf Scotti |42|). 

As the Malliavin calculus has been extended to the case of Poisson measures and SDE's with 
jumps, either dealing with local operators acting on the size of the jumps (Bichteler-Gravereaux- 
Jacod |3] Ma-Rockner|30| Leandre |26] |27| etc.) or based on the Fock space representation of 
the Poisson space and finite difference operators (Nualart-Vives |33] Picard |34] Ishikawa-Kunita 
)20| etc.), it is quite natural to attempt extending the Dirichlet forms arguments to such cases. 
This has been done first by Coquio |16] when the state space is Euclidean then by Denis [18] by 
a time perturbation, see also related works of Privault |36| . Albeverio-Kondratiev-Rockner [T], 
Ma-R6ckner [30]. 

We shall give a general presentation of Dirichlet forms methods for the Poisson measures 
in the spirit of the first approach (Bichteler-Gravereaux-Jacod [3]) which gives rise to a very 
similar situation like in Malliavin Calculus : a symmetric semi-group on the Poisson space and 
a local gradient satisfying the chain rule. With respect to preceding works in this direction we 
introduce a major simplification due to a new tool the lent particle formula |12] which gives the 
gradient on the Poisson space by a closed formula. Thanks to this representation we obtained 
with Laurent Denis several results of existence of density |12] [13] and the method extends to 
C°° results (forthcoming paper). In this lecture I present the method and the main applications 
obtained up to now and I expose new results about the regularity of multiple Poisson integrals 
in connection with the Fock space representation that the Poisson space provides. It is organised 
as follows : 

The functional analytic reasoning. 

Dirichlet forms and non-Gaussian Malliavin calculus — Poisson random measures — Dirichlet 
form on the Poisson space : the lent particle formula. 
Practice of the method. 

Other examples — Applications to SDE's — A useful theorem of Paul Levy. 
Regularity results for multiple Poisson integrals. 

Random Poisson measure and Fock space — Decomposition of D in chaos — Density for 
(h(g),...,I n (g® n )) — Density for (I ni (/f ni ), . . . Jn p UT P )) — Other functionals of Poisson 
integrals — Density of I n (f). 
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2 The functional analytic reasoning. 



Let us first introduce the fundamental notions of the theory of local Dirichlet forms. 
2.1 Dirichlet forms and non-Gaussian Malliavin calculus. 

Let (X, X, v, d, 7 ) be a local symmetric Dirichlet structure which admits a "carre du champ" 
operator. This means that (X, X, v) is a measured space, v is a cr-finite positive measure and the 
bilinear form e[f, g] = \ J j[f, g] du is a local Dirichlet form with domain d C L 2 {u) and carre du 
champ 7 (cf Fukushima-Oshima-Takeda |19| in the finite dimensional case and Bouleau-Hirsch 
in a general setting). The form e is closed in L 2 (y) and the bilinear operator 7 satisfies the 
functional calculus of class C 1 n Lip: 

V/,^d", VF, G of class C 1 D Lip on R n 1 [FU),G(g)]=J2 d i F if)9 j G(9Hfi,9j}- 

ij 

We write always 7 [/] for 7[/, /] and e[f] for e[/,/]. 

The space d equipped with the norm (||.||^ 2 , ^ +e[., is a Hilbert space that we will suppose 
separable. It is then possible to generate the quadratic differential computations with 7 by an 
ordinary differential calculus thanks to the fact that a gradient exists (see Bouleau-Hirsch 
ex. 5. 9 p. 242): there exist a separable Hilbert space H and a continuous linear map D from d 
into L 2 (X, u; H) such that 

• Vued, ||D[ u ]||2 r = 7 [4 

• If F : R —¥ R is Lipschitz then Vn € d, D[F o u] = (F 1 o u)Du, where F' is the Lebesgue 
almost everywhere defined derivative of F. 

• If F is C 1 (continuously differentiable) and Lipschitz from M. d into R (with d S N) then 

d 

Vn = - G d d , fl[Fo U ] = ^T(diFou)D[ui]. 

8=1 

In [11] Chap VII we used for H a copy of the space L 2 (v), but a wide choice is possible depending 
on convenience. 

This differential calculus gives rise to integration by parts formulae as in classical Malliavin 
calculus. For all u € d and v E T> (a) domain of the generator a associated with the Dirichlet 
structure, we have 

- / 7 [it, v\dv = — \ ua[v\du. (1) 



2 . 

The space d n L°° may be shown to be an algebra, hence if iti , 112 € d n L°° 

i J U2j[ui,v]dv = — J uiU2a[v]dv — ^ J ury[u2,v]dv (2) 
Introducing now the adjoint operator 5 of the gradient D, the equality with u 6 d, C7 6 dom <5 

u^C/dr/ = / (_D[u], (3) 

provides for ip Lipschitz 

' ip'(u)(D[u},U) H dv= I tp(u)5Udv. (4) 
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See [7] Chap V to VIII and |9] for applications of such formulae. 

But the Dirichlet structures do possess pecular features allowing to show existence of density 
without using integration by parts arguments. This is based on the following important energy 
image density property or (EID): 

For each positive integer d, we denote by 23(]R rf ) the Borel a- field on M. d and by X d the 
Lebesgue measure on (M. d , B(M. d )). For / measurable f*v denotes the image of the measure v by 
/■ 

The Dirichlet structure (X, X, v, d, 7) is said to satisfy (EID) if for any d and for any M. d -valued 
function U whose components are in the domain of the form 

U4(detj[U,U t ])-u]^:X d 

where det denotes the determinant. 

This property is true for any local Dirichlet structure with carre du champ when d = 1 (cf 
Bouleau |5] Thm 5 and Corol 6). It has been conjectured in 1986 (Bouleau-Hirsch |10| p251) 
that (EID) were true for any local Dirichlet structure with carre du champ. This has been shown 
for the Wiener space equipped with the Ornstein-Uhlenbeck form and for some other structures 
by Bouleau-Hirsch (cf Chap. II §5 and Chap. V example 2.2.4) but this conjecture being at 
present neither refuted nor proved in full generality, it has to be established in every particular 
setting. For the Poisson space it has been proved by A. Coquio )16| when the intensity measure 
is the Lebesgue measure on an open set and we obtained with Laurent Denis a rather general 
condition ( |12| Section 2 Thm 2 and Section 4) based on a criterion of Albeverio and Rockner 
[2] and an argument of Song |44| . The new regularity results that are presented here are based 
on the (EID) property. 

Let us first explain the framework of Poisson measures and the notation of the configuration 
space. 

2.2 Poisson random measures. 

We are given (X, X, v) a measured space. We call it the bottom space. We assume that v is a- 
finite, that for all {x} belongs to X and that v is continuous or diffuse (u({x}) = Vx). 

We consider a random Poisson measure N on (X, X) with intensity measure v. Such a 
random measure is characterized by the fact that for A G X the random variable N(A) follows a 
Poisson law with parameter v(A) and N(A\), . . . ,N(A n ) are independent for disjoint A{. Such 
an object may be constructed on the space of countable sums of Dirac masses on (X, X) (the 
configuration space), by considering first the case where v is bounded where the construction 
is explicit and then proceding by product along a partition of (X, X) (see e.g. |6] or [7] Chap 
VI §3). We denote by (fl,A, P) the configuration space where N is defined, A is the a- field 
generated by iV and P its law. The space (Q, A, P) is called the upper space. 

The following density lemma (cf |12| ) is the key of several proofs. 

Lemma 1. For p G [l,oo[, the set {e~ N ^ : f ^ 0,/ G L 1 (v) D L°° (v)} is total in LP(tl,A,F) 
and{e lN ^: f £ L l (v) D L°°(v)} is total in L P (Q, A, P; C). 

We set N = N - v, then the identity E[(N(f)) 2 ] = f f du, for / G L x (v) n L 2 (v) can be 
extended uniquely to / £ L 2 {y) and this permits to define N(f) for / G L 2 (v). The Laplace 
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characteristic functional is the basis of all subsequent formulae 

E[e mf) ] = e -ni-e*f +l f )d u /GL 2 (z/) (5) 

The creation and annihilation operators e + and e~ well-known in quantum mechanics (see Meyer 
)31| . Nualart-Vives |33| . Picard |34| etc.) will play a central role for calculus on the configuration 
space, they are defined in the following way: 



Vx,w £ Q, = wl {x£suppw} + (w + £ x )l{xtsuppw} 

\/x,w € f2, £ x (w) = wl{ x £ suppw y + (w — £x)l{xesuppw}- 



(6) 



One can verify that for all w G CI, 

£ x (w) = w and £ x {w) = w ~ £ x for iV^-almost all x (7) 

and 

£ x (w) = w + e x and £ x (w) = w for z^-almost all x (8) 
We extend these operators to the functionals by setting: 

£ + H(w,x) = H(e x w,x) and e~H(w,x) = H(e~w,x). 

This extension recommands to be careful with the order of composition since we have for instance 

(£~£ + H)(x, u) = H(x, £+£ x lj) (= H(x, elToj) = e+H) (9) 

It is important to emphasize that since v is continuous the two measures P x v and P^r = 
¥(du)N(uj)(dx) defined on the same sapce (fix X,Ax X) are mutually singular. Computation 
needs to be careful with respect to negligible sets. The next lemma shows that the image of P X v 
by e + is nothing but Pyv whose image by e~ is P x v : 

Lemma 2. Let H be A® X -measurable and non negative, then 

E / e + Hdv = E / HdN and E / e~HdN = E / Hdu. 



We will encounter also another notion, sometimes called a "marked" Poisson measure associ- 
ated with N, which needs here a rigorous construction. 

We are still considering N the random Poisson measure on (X, X, v) and we are given an 
auxiliary probability space (R,lZ,p). We construct a random Poisson measure iV p on (X x 
R, X ® H, v x p) such that if iV = ^ i e Xi then N Q p = ^ £ (xi.ri) where (/j) is a sequence of 
i.i.d. random variables independent of N whose common law is p. 

The construction of N p follows line by line the one of N. Let us recall it. We first study the 
case where v is bounded and we consider the probability space (N, P(N), P u (x)) x %i > 
where P u {x) denotes the Poisson law with parameter v(X) and we put 

Y 
N = e x . , (with the convention = 0) 
i=i l 

where Y, x±, ■ ■ ■ , x n , ■ ■ ■ denote the coordinates maps. We introduce the probability space 

(£l,A,P) = (R,K,pf\ 
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whose coordinates are denoted by n,--- , r re , • ■ • . On the probability space (N,"P(N), P v [x)) x 
(X, X, ^nc) ) N * x ^> we define the random measure N Q p = ^2i=l £ (xi,n)- ^ ^ s a P°i sson 
random measure on X x R with intensity measure u x p. For / 6 L 1 ^ x p) 

E[[ fdNQp]= [ ([ f{x,r)dp{r))N(dx) P - a.e. (10) 

JXxR JX JR 

and if / G I?(v x p) 

E[( / /cZiV p) 2 ] = ( / / fdpdN) 2 - I (I fdpfdN +11 fdpdN, (11) 
JxxR Jxjr Jx Jr Jxjr 

where E stands for the expectation under the probability P. 

If v is fj-finite, this construction is extended by a standard product argument. Eventually in 
all cases, we have constructed N on (fl, A, P) and N Q p on (Q,A,¥) x (fi, .A, P), it is a random 
Poisson measure on X x ii with intensity measure v x p, and identities fjlOf) and HI 1 [) generalize 
as follows: 

Proposition 3. Let F be an A (g> X <S> 1Z measurable function such that E f XxR F 2 dudp and 
E ffjifx \F\dv) 2 dp are both finite then the following relation holds 

E[( / FdN p) 2 } = ([ [ FdpdN) 2 - [ ( / FdpfdN + [ [ F 2 dpdN, (12) 

JXxfi JXJR JX JR JXJR 

in particular if F is such that J Fdp = P x v-a.e., then E[(j XxR FdN p) 2 } = j x J R F 2 dpdN. 

Proof. Approximating first F by a sequence of elementary functions and then introducing a 
partition (Bk) of subsets of X of finite v- measure, this identity is seen to be a consequence of 
(HU). □ 



Let us take the opportunity to state two formulae that we didn't mention in our preceding 
articles, and which may be useful in some context. Let F be measurable as in Prop [3] and say 
< F < 1 then 

E exp J log F dN p = exp J (log J F dp) dN (13) 



E J FdN p= I (J Fdp)dN (14) 
whose proofs follow the same lines as the construction of N Q p and Prop [3] 

2.3 Dirichlet form on the Poisson space : the lent particle formula. 

Now, after these notions related to the pure probabilistic Poisson space, we shall assume we 
have on the bottom space a Dirichlet structure {X, X , v, d, 7) as defined in section 2.1. And we 
attempt to lift up this structure to the Poisson space in a natural manner. This may be done in 
several ways (see e.g. the introduction of |12]). The method we will follow is not the simplest, 
we choose it because it enlightens the role of operators e + and e~ in the upper gradient. 
First, thanks to (|5|) we obtain the following relation: for all / £ d and all h G T>(a), 



E 



i *V ) [N(a[h]) + l i N(>y\f,h]) 



0. (15) 
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This relation and the explicit construction which may be done when v is a bounded measure (cf 
[6]) suggest a candidate for the generator of the upper structure. 
Let us consider the space of test functions 

V = £{e^ (/) with / G V(a) D L X (u) et 7 [/] G L 2 (v)}. 

and for U = Y, p \e ifi ^ p) in V , let us put 

Ao[U] = y £^ i * iU) (iN(a[f p }) - |JV( T [/ P ])). (16) 
p 

The procedure to show that Aq is uniquely defined and is the generator of a Dirichlet form 
satisfying the hoped properties, has two steps : first to construct an explicit gradient, then to 
use Friedrichs' property. 

Gradients. 

We will suppose as in section 2.1 that the bottom structure possesses a gradient that we denote 
from now on (•) . For convenience we assume it satisfies the following properties 

• constants belong to d/ oc (see Bouleau-Hirsch Chap. I Definition 7.1.3.) 

1 G d loc which implies 7[1] = and l b = 0. (17) 

• (.)'' is with values in the orthogonal subspace Lq(R, 1Z, p) of 1 in the space L 2 (R,7Z, p). 
This condition is costless since for the gradient only the Hilbert structure of H matters. From 
now on we denote this gradient (.)'". 

We take for candidate of the upper-gradient for FG^o the pre-gradient 

F s = J e~((e + F) b ) dN p. 

where N p is the Poisson measure N "marked" by p as defined in section 2.2. 
Let us remark that thanks to Prop 3 and (j 1 7[) we have 

E[(/ e-((e + Ff)dN Q p) 2 ]= [ e-(-f(e + F))dN F-a.e. (18) 

JXxR J 
For / G V{a) n L l (m), j[f] G L 2 , we have e ifi{f) G V and 

(e^(/))» = y e^iiffdNQp 

what yields on Dq: 

E[F»G»] = ]T \ p JT q e^f^N( 7 (f p , g q )) (19) 
P,q 
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Friedrichs' argument. 

This enables us to show that the representation (11 6|) does not depend on the expression of U 
and that Aq is indeed a symmetric negative operator on the dense subspace T>q of L 2 (P) so that 
Friedrichs' argument applies (see [TT] p. 4 or [7] Lemma III. 28 p. 48) : it can be extended to a self 
adjoint operator which may be proved to generate a Dirichlet form with domain B admitting a 
carre du champ T with a gradient extending (.)". 

It remains only a technical point to verify: the fact that T>q be dense in L 2 (P). This is not 
obvious because of the condition j[f] G L 2 (v) that we need in T>q in order Aq take its values 
in L 2 (P). In |12| we called it bottom core hypothesis (BC), it is not a real constraint in the 
applications. We can state (cf |12| ) : 



Theorem 4. The formula 



VFGl, F*= / e~((e + Fy)dN © p, 

JExR 



'ExR 

extends from T>q to B, it is justified by the following decomposition : 

FeB £ t^ 1 e + F - F gD^' £-((e+F) b ) G L 2 (F N x p) d( ^ p) F» G F 2 (P x P) 

where each operator is continuous on the range of the preceding one and where Lq(Ptv X p) is 
the closed set of elements G in L 2 (Pjv x p) such that J R Gdp = P7v-a.s. Furthermore for all 
F G 



F[F] = E^) 2 = [ e~ 7 [e + F] dN. 
Je 



This main result — that we call the lent particle formula — implies the validity of a functional 
calculus for the obtained Dirichlet structure (f2, A, P, B, T) on the Poisson space that may be 
sketched as follows: 

Let be H = . . . ,F n ) with $ G C 1 n Lip(M n ) and F = (Fi, . . . , F n ) with Fj G B, we 

have : 

a) l[e + H] = Zij me+F^ie+F^ie+Fi^+Fj] P x ^-a.e. 

b) e- 7 [e+H] =^2 ij & i (F)& j (F)e- 1 [e+F i ,e+F j ] P^-a.e. 

c) T[H] = Je--f[e + H}dN = $' i (F)^.(F)/e- 7 [£+F i ,e+F ? -]cZiV P-a.e. 
Remark 5. Let F G B, by the theorem applying formula (jl 3[) to F" gives 

EexpF tt = Eexp j e~(e + F) b N Q p = exp ^ ^log ^ expe"(e + F) b (ip^ diV 

= exp J (e~ log y exp (e+F^dp) dN (20) 

what may yield the characteristic function of the law ofF" under PxP: if we put J exp(me" 1 " F)^ dp 
exp $ (it) we obtain 

KEe iuFt =Eexp J e~^(u)dN. 
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Example 1. 

Let F( be a centered Levy process with Levy measure a integrating x 2 and such that a local 
Dirichlet structure may be constructed on M\{0} with carre du champ -y[f] = x 2 f' 2 (x). With 
our notation (X, X,v) = (M+ x ]R\{0}, Borelian sets, (it x a). 

We define the gradient b associated with 7 by choosing £ on the auxiliary space (Cl,A, P) 
such that J ^(r)dr = and L £ 2 (r)dr = 1 and putting = xf'(x)^(r). 

The operator b acts as a derivation with the chain rule (y(/)) = f'{f)-f^ (f° r E C 1 H Lip 
or even only Lipschitz). 

N is the Poisson random measure associated with Y with intensity dt x a such that Jq h(s) dY s - 
J 1[ j](s)h(s)xN(dsdx) for h G L 2 oc (R+). (These hypotheses imply 1 + AY S / a.s.) 

Let us study the existence of density for the pair (Yj, £xp(Y)t) where £xp(Y) is the Doleans 
exponential of Y. 



£xp(Y) t = e Yt H(l + AY s )e~ AY \ 

s ^ t 

1°/ We add a particle (a, y) i.e. a jump to Y at time a ^ t with size y : 

e+ ay) (£xp(Y) t ) = e Yt+ y J] (1 + AY s ) e - Ay *(l + y)e^ = £xp(Y) t (l + y). 



s < t 



2°/ We compute 7[e+£ xp(Y) t ](y) = {£xp(Y) t ) 2 y 2 . 
3°/ We take back the particle : 

e-jle+SxpWt] = (£xp{Y) t (\ + y)- x ) 2 y 2 
we integrate in TV and that gives the upper carre du champ operator (lent particle formula): 

F[£xp(Y) t ] = Jj 0>t]xM (Sxp(Y) t (l + yr 1 ) 2 y 2 N(dady) 
= E a {£MYW + A^)" 1 ) 2 AY 2 . 

By a similar computation the matrix T of the pair (Y~t, £xp(Yt)) is given by 
r = ( 1 £xp(Y) t (l + AY a )^ \ 

= ^ t V £x P (Y) t (i + aYq,)" 1 (£xp(y) t (i + aYq,) -1 ) 2 y a • 

Hence under hypotheses implying (EID) the density of the pair (Yt,£xp(Yt)) is yielded by the 
condition 

dhaC ({£xp(Y) t (l + AY a )^) ^ JT )=2 

where JT denotes the jump times of Y between and t. 
Making this in details we obtain 

Let Y be a real Levy process with infinite Levy measure with density dominating a positive 
continuous function ^ near 0, then the pair (Y t ,£xp(Y) t ) possesses a density on R 2 . 
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Example 2. 

Let Y be a real Levy process as in the preceding example. 

Let us consider a real cadlag process K independent of Y and put H s = Y s + K s . Putting 
M — sup s < j. H s and computing successively (e^ikf), 'yfe^iW] and applying the lent particle 
formula gives 

Proposition 6. // a(IR\{0}) = +00 and if P[sup s < t H s = Hq] = 0, the random variable 
sup s < t H s possesses a density. 

It follows that any real Levy process X starting at zero and immediately entering whose 
Levy measure dominates a measure a satisfying Hamza's condition ([19] pl05) and infinite, is 
such that sup s < t X s has a density. 

Example 3. Levy's stochastic area. 

This example will show that the method can detect densities even when both the Malliavin 
matrix is non invertible and the Levy measure is singular. 

Let X(t) = (Xi(t), X2(t)) be a Levy process with values in M. 2 with Levy measure a. We 
suppose that the hypotheses of the method are fulfilled, we shall explicit this later on. 

Let us consider first a general gradient on the bottom space : 

f = /& + fk2 

where f[ = J||-, and £1, £2 are functions defined on M. 2 x R which satisfy: J R £i(-, r)p{dr) = 
J R £,2{-,r)p{dr) = 0, f R t%(x 1 ,x 2 ,r)p(dr) = a 11 (x 1 ,x 2 ), f R £i(x 1 ,x 2 ,r)£ 2 (x 1 ,x 2 ,r)p(dr) = a 12 {xi,x 2 ), 
fn£,2( x i' x 2)p(d r ) = 022(^1,3:2), so that 

7[/] =«n/i 2 + 2012/1/2 + a 22 /2 2 - 
Let us consider the following vector involving Levy's stochastic area 

V = (X 1 (t),X 2 (t), f X 1 {s-)dX 2 {s) - f X 2 {s-)dX x {s)). 
Jo Jo 

We have for < a < t and x = {x\,x 2 ) € M 2 , 

£ t,x) V = V + {xi,x 2 ,X 1 (a„)x2 + x x {X 2 {t) - X 2 {a)) - X 2 (o_)x 1 - x 2 {X x {t) - X x {a)) 

= V + (x 1 ,x 2 ,x 1 (X 2 (t) - 2X2(0)) - x 2 (X 1 (t) - 2Xi(o))) 
because e + V is defined P x vxda-a.e. and v x da is diffuse, so 

(e+vf = (6,£ 2 ,£i(x 2 (t) - 2X2(0)) - 6(Xi(i) - 2X1(0))) 

and 

(on 012 Aaii — Bai 2 \ 

012 o 2 2 ^012 - -Bo 2 2 

Aaii ~ Bai2 Aa X2 - Ba 2 2 A 2 au - 2ABa X2 + -B 2 o 2 2 / 

denoting A = (X 2 (i) - 2X 2 (o)) and B = (Xi(t) - 2Xi(o)). 
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This yields 



e~ A = X 2 (t) - AX 2 (a) 



2X 2 (a_) 



let us denote it A 



e~B = X 1 (t) - AXi(a) 



2Xi (a_ ) 



let us denote it B 



and eventually 



m = J2 



( 



an(AX a ) a 12 (AX a ) 
~ a 2 2(AI Q ) 



A 2 



A ail {AX a ) - Ba 12 (AX a ) 
Aa 12 (AX a ) - Ba 22 (AX a ) 
au(AX a ) - 2ABa 12 (AX a ) + B 2 < 



a 22 (AX a ) 



) 



the symbol ~ denoting the symmetry of the matrix. 

Considering the case a± 2 = let us take the Levy measure of (X±,X 2 ) expressed in polar 
coordinates as 



with g locally bounded and such that it dominates a continuous and positive function near 0. 



Then V = (X^t), X 2 (t), f* Xi(s-)dX 2 (s) - f* X 2 (s-)dXx(s)) has a density (and condition (0.4) 
of |14j or of )34j propl.l are not fulfilled). 



The Levy measure of (Xi, [Xi]) is carried by the curve x 2 = x\. We have \(xi,x 2 ) = 2x\. 
We arrive to the sufficient condition : V has a density as soon as the Levy measure of X\ is 
infinite and satisfies hypotheses for (BC) and (EID). (cf [12] and |13|). 

3 Practice of the method. 

Computation with the lent particle formula. 

The presence of operators e + and e~ in the lent particle formula (Thm [4]) which exchange the 
mutually singular measures Pat and Fxv, requires to be more careful than in the usual stochastic 
calculus where all is defined P-a.s. We make some remarks and give some examples to help the 
reader to become familiar with this tool. 

The lent particle formula extends to ID>i oc . 

The space D; oc is a remarkable specific feature of local Dirichlet forms with carre du champ : 
the carre du champ operator extends to functions locally - in a measurable sense - in D (cf |11] 
Chap I §7.1). 

We denote D; oc the set of applications F : £1 1— > 1R such that there exists a sequence £l n G A such 
that U n Q n = £1 and 3F n E D with F = F n on tt n . 

The fact that (EID) is always true for d = 1 (cf [5]) shows that, for F € ©z oc , T[F] is uniquely 
defined and may be evaluated by r[_F n ] on f2„. The operator jj extends to B; oc by putting F^ = f| 
on £l n . For F in 0; oc , the formulae 



v(dp,M)=g(9)de.l ]0tl[ (p) 



dp 



P 





resume a computation done on each 11 
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Negligible sets. 



As it was recalled above at the beginning of section 3, it is recommended to write down the 
negligible sets at each equality e.g. 

e+(Nf) = Nf + f Fxu-a.e. 

e~(Nf) = Nf-f F N -a.e. 

e + (e iNf g) = e iNf e if g F x v-a.e. 

£~{e lRf g) = e iRf e~ if g F N -a.e. 

Remark 7. Let us observe that if H{u,x) = G(co)g(x) where G is defined P-a.s. and g v-a.e. 
then H belongs necessarily to a single class P^v-a.e. So that we may apply to H both operators 
e + and e~ without ambiguity. This will be used further about multiple Poisson integrals. 



A simplified sufficient condition. 

Theorem H] gives a method for obtaining T[F] for F 6 D or F E B n , then with the hypotheses 
giving (EID) it suffices to prove det T[F] > P-a.s. to assert that F has a density on W 1 . Let 
us mention a stronger condition which may be also useful in some applications. By the following 
lemma that we leave to the reader 

Lemma 8. Let M a be random symmetric positive matrices and fi(da) a random positive mea- 
sure. Then {det / M a fi(da) = 0} C {/ det M a (i(da) = 0}, 

it is enough to have J det e~( / y[e + F])dN > P-a.s. hence enough that det e~ ( , y[e + F]) be > 
Pjv-a.e. We obtain, by lemma 2, that a sufficient condition for the density of F is det ^y[e + F] > 
P x v x dt-a.e. (or equivalently that the components of the vector (e + F)^ be P x v x dt-a.e. 
linearly independent in L 2 (p) ). 



The energy image density property 

(EID). We gave in Bouleau-Denis |12] general conditions on the bottom structure (X, X, u, d, 7) 
to satisfy (EID) and for this property to be lifted up to the upper space (O, A, P, D, T). Here are 
these conditions in a simplified form: 

Proposition 9. Suppose (X, X, v) = (R , B(M, d ), k{x)dx) with k continuous on an open set of 
full Lebesgue measure and suppose the carre du champ operator is defined on the test functions 
C|? infinitely differentiable with compact support by the formula 

£&(*W(aO07/(aO (21) 

where £ is locally bounded and locally elliptic i.e. for every compact K there are constants 
Ck < 00 and ck > such that Vx G K, Vc £ M. d Ck\c\ 2 ^ Yli j=i dj( x ) c i c j ^ ck\c\ 2 , then the 
bilinear form 

e[u, v] = - / £ij(x)dju(x)djv(x)k(x) dx. (22) 
2 Jw _ 

defined on is closable and its closure defines a Dirichlet form (e, d) with carre du champ 
given by (|2ip . and this structure satisfies (EID) and (BC). 
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It is useful for many examples to remark that the preceding case allows to extend (EID) and 
(BC) to situations where v is singular w.r. to Lebesgue measure. 

Let (RP\{0},£(IR p \{0}),^d,7) be a Dirichlet structure on W\{0} satisfying (EID). Let 
U : W\{0} ^ R q \{0} be an injective map (p < q) such that U G d q . Then U*v is cr-finite. If 
we put 

du = {<p G L 2 (U*v) : ip o U G d} 
eu[ip] = e[<poU] 

r I dU* hlfoU}. v) 

yum = — du t v 

then the term (M 9 \{0}, #(M 9 \{0}), U*v, djj, <yu) is a Dirichlet structure satisfying (EID). Addi- 
tional regularity assumptions make U transport also property (BC). 

Now it is possible to lift up (EID) from the bottom to the upper space if two conditions are 
fulfilled. First to be able to share the bottom space on a partition of sets of finite z/-measure. 
Second that the obtained Dirichlet structures are such that any finite product satisfies (EID). The 
precise formulation is given in Bouleau-Denis |12| Section 4. This covers all cases encountered 
in practice. 



3.1 Other examples. 

Example 4. Nearest point of the origin. 

This example shows the quickness of the method which has, in some sense, to be paid by the 
care to put on negligible sets. 

Let us take for the bottom space (R d , B(R d ),u, d, 7) satisfying (BC), assuming the identity 
map j on W d belong to d d and 7[|j|] > 0, the measure v being infinite, possibly carried by a 
surface or a curve. Let us consider the functional H defined on (fl, A, P) 

H(oj) = inf Ixl. 
xesupp(w) 

The inf is reached because the measure v is a- finite. We have 

e+H = \x\ A H P x zA-a.e. 

We will suppose that the measure v does not charge the level surfaces of \x\ i.e. the spheres 
centered at O. Then for fixed us, x i— > £%H belongs to d and we have 

(etHf = (|j|)%| ^ H = (\j\)\\<H Fxn p-a.e. 

The two functionals lu-i ^ u and hj\<H equal P x v-a.e. do have the same image by e~ P^v-a.e. 

1 \x\ < H(e^cu) = l \x\<H(e*u) ^N-a.e. 
and the lent particle formula gives 

r[H] = J y[\j\}(^\ x \ <H(£ - u) N(u>,dx) = J l[\j\](x)l\ x \ <H{e - u) N(u,,dx) P-a.s. 

Now this integral is easily seen to be equal to 7[|i|](a?o( w )) where xq is the P-a.s-unique point 
achieving the minimum of the distance of the support of u to the origin. Thus we obtain the 
quite natural result that as soon as v doesn't charge the spheres, H possesses a density. 

As in several other examples, the result could be extended to the case where is(M. d ) be finite 
by conditioning by the event {A^B^) ^ 1}. 
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Example 5. Gas of Brownian particles. 

This is an extension of the preceding example to infinite dimensional setting. We consider a gas 
of Brownian particles in M 3 . Each particle is independent, the initial positions are distributed in 
R 3 along a Poisson measure with uniform intensity. We study the lowest distance of a particle 
to the origin during the time interval [0, 1]. 

A) Let us begin with some properties of extrema on the Wiener space. Let be given a Brownian 
motion Bt = (Bj,Bf,Bf) starting at zero, the Wiener space being endowed with the Ornstein- 
Uhlenbeck structure. We adopt — only in this paragraph A) — the following notation for this 
structure (W, W, m, B, T) and we use the gradient with values in L 2 (W, W, rh) where (W, VV, to) 
is a copy of (W, W, to) defined by 

(f f{t)-dB t f= f f(t)-dB t V/ = (/ 1 ,/ 2 ,/ 3 )GL 2 ([0,l]). 
Jo Jo 

If x G M 3 is fixed and ^ 0, the random variate 

K(w) = inf \x + BA 
te[o,i] 

is strictly positive and in D, by the argument developped by Nualart-Vives |32| . using the fact 
that the set of Brownian paths which reach several times the minimum is negligible, we obtain 

jAl M (*i + B^ w) (w))B^ (w) (vj) + (*2 + B 2 T[w) {w))B 2 T{w) {w) + (x 3 + B* {w) (w))B* {w) (w) 

K»(W,W) = : — 

\x + B T{w) {w)\ 

where T(w) = mi{t G [0, 1] : \x + B t {w)\ = K(w)}. 
It follows that 

T[K] = E[(K$) 2 } = T > a.s. ifx^O. 

B) Let us come back to our usual notation. For the bottom space we take (X, X, v) = (M 3 x 
W, jB(M 3 ) x W, A 3 x m) where A 3 is the 3-dimensional Lebesgue measure, that we equip with the 
product Dirichlet structure of the zero form on M 3 and the O-U-form on the Wiener space. The 
structure (X,X,v,d,j) is thus naturally endowed with a gradient induced by the gradient used 
in part A) and that we denote now b as usual, it is with values in L 2 (m). The hypothesis (BC) 
is fulfilled. 

We construct the upper structure (fi, A, P, D, T) which describes a gas of Brownian particles. 
We denote (x, w) the current point of X and we consider the functional 

H(u)= inf \x + B t (w)\. 

t G [0, 1] 
(x, w) G supp to 

We apply the lent particle method : 

s+ w) H = ( inf \x + B t (w)\) AH 
yx,w) te[o,i] 

Here the measure A 3 x m does not charge the level sets of (inf tg [o,i] \ x + Bt(w)\) and we have 

(e+H'f = (inf tg[0jl] \x + £ t (u>)|) b l {(infig[M \ X+Bt { w )\) ^ H} 

= ( mf tG[o,i] \x + B t (w)\)h {iin{te[0A] \ x+Bt (w)\)<H} P x v x m-a.e. 
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what gives putting rj(x,w) = inf fe [ .i] \ x + Bt(w)\ and a = (x,w) 

T[H] = J l[v}\ v{a) ^H { s-.)}N^,da) 

and this is equal to 7(77] taken on the unique Brownian particle which yields the minimum. 

Since by the part A) above this quantity is strictly positive, we can conclude that H possesses 
a density. 

The argument extends to the case where the point taken as origin is itself moving determin- 
istically or as an independent diffusion process. 

Example 6. Integral of a Levy process. 

let Yt be a Levy process with values in M. d with our usual hypotheses that the Levy measure a 
carries a Dirichlet form such that hold (BC) and (EID). Let us suppose in addition for simplicity 
that a integrates |x| 2 , that Y is centered without Brownian part and that the coordinate maps 
Xi are in d with 7^,^] = XiXjSij. 

Let be g € C 1 n Lip from M. d into itself and let us consider the d-dimensional functional 

H= f g(Y t )dt 
Jo 

which writes also H = Jq g{f l[o,t] (s)y N(dyds))dt if N denotes the Poisson measure associated 
with (Y). For < a < 1 

e t,y) H = Io9(Yt)dt + j l a g{Y s + y)ds P x „-a.e. 
(s^ y) Hf = j 1 a Dg(Y s +y)d S -j\y) P x v x p-a.e. 

where Dg is the Jacobian matrix of g and j the identity map on M. d . Then 

e-(e+ jS/) tf) b = jlDg(Y s )ds-j\y) F N x p- a .e. 

"m] = E a *i J« D 9(Y 8 )d8 7\j,3*](AY a )£D*g(Y 8 )d8 P-a.s. 

the matrix r y\j,j t ](AY a ) is the d x (i-matrix whose diagonal is composed of the squares of the 
jumps ((AY^) 2 , . . . , (AY^) 2 ). If the images by the coordinate mappings of the Levy measure are 
infinite, and if the Jacobian matrix Dg is regular, then H has a density on M. d . See |25| and [3] 
for related results. 

Example 7. Generalized Ornstein-Uhlenbeck processes. 

Let (£,77) be a 2-dimensional Levy process starting from (0,0). The process 

X t = e^{x+ f e-^dr] s ) t^O xE« 
Jo 

is a homogeneous Markov process called generalized O-U process driven by (£,77) (cf |15|). It is 
possible to see by the classical Malliavin calculus that if (£, 77) possesses a Brownian part then Xt 
has a density. We exclude this case now and suppose that the Levy measure carries a Dirichlet 
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form satisfying (BC) in order to apply the method (without care of (EID) because Xt is one 
dimensional) . 

Let us begin by computing T[Xt] by the lent particle method. 
Let (a, £,rj) denote the current point of X = R + xlxl, 



(e + X t f 
e-{e + X t f 



X + J[0,t] e 



6- = 6- + ei 



d(rj s + r/l Q <; . 
£ b + /[o,a] e-«- + f 1 «<')d77 s ^ + e-*- (?? b + ??£ b 



+ I[o,a] e ~ is ' d Vst + e-^-rf 



Let j be the identity map on M?, so that j* = (£' ,rf) we obtain 
*t 



T\X t 







+ 

now putting V(lj, x, a, £, 77) 



[0,q] 



drjs 



x + / a e f 



drj s 



N(dad£dr]) 



dr] s 



this writes 



T[X t ] = e 2 ^ Yl VHjJ'KMa^^V 



the sum being taken on the jump times of the process ((,t,Vt)- Starting from this relation we 
discuss several cases : 

1) First case detj[j, j*] > 0. 

Since (EID) does not matter Xt being real valued, the only condition is that the Levy measure 
a of (€t)Vt) be infinite and carry a local Dirichlet structure (R 2 \{(0, 0)}, £>(R 2 \{(0, 0)}), a, d, 7) 
satisfying (BC) and such that j £ d/ oc and det > cr-a.e. 

No necessary and sufficient condition is known for this which would extend the Hamza condition 
to dimension 2, but we see by Prop 7 that this will be fulfilled as soon as a has a continuous 
density. 

2) The case where ^ and rjt are independent. 



The measure a is carried by the coordinate axes, 7 [7, j*](a, b) 



ip{a)l b=0 
ip(b)l a=0 



and 



E 

a sS t 



+ / e-^-di h ) 2 V (AC a )l Ava=0 + e- 2 ^-^(A Va )l Aia 
J[0,a] 



If the Levy measure of (rjt) is infinite and if ip > then Xt has a density. 

If the Levy measure of (rjt) is finite, then if (p > and if the Levy measure of (£t) is infinite Xt 
has a density as soon as x + Jj a i e~^ s ~di] s does not vanish for little a hence as soon as 1 / 0. 
In this case of independence it is also possible to use the representation in law (cf |15] Thm 

3.1) 



e^x + 



drjs 



16 



applying the lent particle method to the right hand side. That gives a little faster the same 
conclusion. 

3) The case where a is carried by a curve. 
We sketch only this case which involves a parametrization. Let M be a Poisson measure on R + 
with CT-finite intensity measure m and let be given a map <E> : u h- > (f(u),g(u)) 6 M 2 \{(0, 0)} 
such that we obtain our Levy process by image : 

= a $*M = 2V 

$ being injective from R + into M 2 \{(0, 0)} and such that lim a _ i . 00 <&(a) = (0, 0). On BL + we start 
with a Dirichlet structure (M+, £?(IR + ), to, d, 7). We assume the identity J € d/ oc and / and g of 
class C 1 fl Lip. 

f' 2 f'sf 



We have 7 [7, j*] = ( ^ , ^ f 2 ) an d we obtain 



r[x t 



1 2 



/'(a)(s+ / e 5 *-d7/ s ) +^(a)e" 



j[J](a)M(dada). 



Let us suppose the Levy measure of (^aVt) infinite, i.e. m infinite, and lim (1 _ j . 00 (/ / (a), g'(a)) 
exist and be equal to (v\,V2) 7^ (0,0). 

T[Xt] = for some oj would imply v\x + 1)2 = what can be realized only for one value of 
x. The reasoning may then be improved by considering the behaviour at the neighborhood of 
another time q.q. 

Example 8. Interaction potential. 

Several forms of interaction potential are encountered in physics for an infinite system of inter- 
acting particles: exp{-/3 m(X t - X 3 )} , ap n Uij g(\Xi ~ Xj\) or exp^- a(^i)«(^)KI^ - 
Xj\)} etc. 

Let us consider the functional $ = J tp(x)ip(y)ip(\x — y\ 2 )N(dx)N(dy) where the functions <p 
and ij) are regular, ip(0) = 0, N being a random Poisson measure on IR 3 . 

After computing as usual e^"$, (e + $) 1 ' and e _ (e + $) 1 ', the lent particle theorem gives 



r[$]= / v{x)Hi,f\v{x) N(dx) 

where j is the identity on R 3 and V(x) is the column vector 

V(x) = j (2<p(a)if>(\x - a] 2 )V<p(x) + 4(p(x)(p(a)i/j' (\x - a\ 2 )(x - a)) N(dx). 
If the bottom structure is such that j*] may be chosen to be the identity matrix, we have 



r[*] = JI F(x,y)N(dy) 



N(dx) (23) 



with F = [2ip(\x - y\ 2 )Vip(x) + A<p{x)tp'{\x - y\ 2 ){x - y)]tp(y). 

In order to study the positivity of r[$], we will use the following lemma (due to Paul Levy 
1931) on which we will come back in the next section. 
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Lemma 10. Let / be measurable on the bottom space such that J \ f\ A \dv < +00. 
If v{f 7^ 0} = +00 then the law of N(f) is continuous. 

That gives us the following result 

Proposition 11. If F is such that (i) 3G € L x (v) : \F(x,y)\ ^ G(y), (ii) Vy 14 F(x,y) is 
continuous, (iii) Vx u{F(x, .)} = +00, then (J" F(x,y) N(dy) 7^ 0) P-a.s. 

Proof. For u) outside a negligible set F(x, .) is bounded in modulus by an integrable function 
for N(u,dy), hence x 1— >■ J F(x,y)N(dy) is continuous by dominated convergence, hence the 
set {x : J F(x, y)N{dy) 7^ 0} is open; by the property (iii) and the lemma this set contains a 
countable dense set, hence all the space. □ 
It follows that if the bottom structure satisfies (BC) $ has a density. 

3.2 Application to SDE's. 

Let d € N*, we consider the following SDE : 

X t = x+ [ [ c(s,X s -,u)N(ds,du) + f a(s,X s -)dZ s (24) 
Jo Jx Jo 

where x € R d , c : M + x R d x X ->• M d and cj : R+ x R d -)• M dxrt , Z is a semi-martingale and iV 
a compensated Poisson measure. 

The lent particle method allows to apply the machinery of Malliavin calculus faster than 
usual and under a set of hypotheses that express the Lipschitz character of the coefficient and 
some other regularity assumptions for the details of which we refer to |13] , 

Let us emphasize that applying the method to SDE's uses reasoning in complete functional 
spaces in which may be computed and solved the stochastic differential equations giving the (j 
of the solution. This takes full advantage of the fact that the lent particle formula is proved not 
only on a set of test functions but on the space B itself. 

In )13| applications are given to McKean-Vlasov type equation driven by a Levy process and 
to stable like processes. 

Example 9. A regular case violating Hormander conditions. 

The following SDE driven by a two dimensional Brownian motion 

' X\ =z 1 + f*dBl 

< Xl =z 2 + J t Q 2XldBl + j t dBl (25) 
k X? =z 3 + tixtdBl + 2f*dBl 

is degenerate and the Hormander conditions are not fulfilled. The generator is A = \{Ul+U$)+V 
and its adjoint A* = \{Uf + E7|) - V with U x = £- + 2xi^ + Xl ^, U 2 = ^ + 2^ and 
V = — \ ~§^l~ The Lie brackets of these vectors vanish and the Lie algebra is of dimension 

2: the diffusion remains on the quadric of equation \x\ — x 2 + \x% — ft = C. 
Let us now consider the same equation driven by a Levy process : 



IS 



Z\ =z 1 + tidY s 1 

Z 2 =z 2 + ti2ZldY s 1 + f*dY s 2 (26) 
Zf =z 3 + f*ZldY s 1 + 2f*dY? 

under hypotheses on the Levy measure such that the bottom space may be equipped with the 
carre du champ operator j[f] = y\f 2 + y 2 f 2 satisfying (BC) and (EID). Applying the lent 
particle method is as usual and shows easily that if the Levy measures of Y 1 and Y 2 are infinite 
Z t has a density on M 3 . See |12j for details. The regularizing property is related to the fact that 
equation (|26|) is not under the canonical form in the sense of Kunita [23] |24| . The next example, 
on the contrary shows a Levy process in M 3 living on a hyperbolic paraboloid. 

Example 10. 

For a£l 3 , let us consider the diffusion solution of 



X t = a+ f U 1 (X s )odB 1 s + [ U 2 (X s )odBl 
Jo Jo 



where B = (B 1 ^ 2 ) is a standard Brownian motion with values in IR 2 , integrals being in the 
Stratonovich sense, and vectors U\ and XJ 2 being given by 



x\x\ - a x 2 x 3 \ / x 2 x\ + a xix 3 

Ui(x) = | x 2 x\ + a xix 3 U 2 (x) = xixl - a x 2 x 3 

x 3 (a 2 +x 2 ) J \ x 3 (a 2 +x 2 ) 



with ao = a 2 + a 2 . — a 2 . 

Then the diffusion (Zt) remains on the quadric of equation 

x\ + x\-x\ = al. (27) 

Now let us consider two independent Levy processes (Y^), (Y 2 ) and the equation 



Z t = a + f U^Z^dY} + [ U 2 {Z s _)dY: 
Jo Jo 



(28) 



the Markov process with jumps Z remains on the hyperbolic paraboloid (|27p as seen by applying 
Ito formula. This is due to the fact that the HP is a ruled manifold and at each point of it 
the jumps of Z are in the direction of either generatrix crossing at this point. Equation f)28|) is 
canonical in Kunita's sense. Using a map from the HP to M 2 the method allows to show the 
density of the law of Zt w.r. to the area measure on the HP. 



3.3 A useful theorem of Paul Levy. 

It is the occasion to rectify a historical injustice about the remarkable article of Paul Levy "Sur 
les series dont les termes sont des variables eventuelles independantes" which appeared in Studia 
Mathematica in 1931 [28J. This article is almost never cited up to now (today the search engins 
do not mention any citation of this article) and the textbooks of K.I. Sato |41| and of J. Bertoin 
do not quote it. One of his theorems, that we recall below, is generally attributed to Hartman 
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and Wintner "On the infinitesimal generator of integral convolutions " Amer. J. Math. 64, (1942) 
273-298, which was published ten years later. 
Paul Levy's results^ may be stated as follows: 

Theorem 12. Let X n be a sequence of independent real random variables such that the series 
^2 X n converges almost surely. 

a) If for any sequence of constants (a n ), ^P{X n ^ a n } diverges, ^2X n has a continuous 

law. 

b) If there is a sequence (a n ) s.t. ^P{X n ^ a n } converges and if the lower bound of the total 
mass of the discrete part of the laws of the X n 's is zero, then the law of ^2 X n is continuous. 

It follows from this theorem that any process with independent increments whose Levy mea- 
sure in infinite has a continuous law. In the framework of random Poisson measures it gives 
easily Lemma [TOl above. 

Remark 13. If / £ and v{f ^ 0} = +oo then the law of N(f) is continuous but 

its characteristic function does not necessarily tend to zero at infinity, in other words is not 
necessarily a Rajchman measure (cf |37] |38] or |8]) . This gives an easy way to construct 
continuous measures which are not Rajchman. Let m = f*v, m is cr-finite and integrates x i— > \x\. 
Since Ke luN f = e-^ 6 "" ~ 1 + lu f) dl/ the law of Nf is Rajchman iff lim^ u ^ +00 J(l — cos ux)m{dx) = 
+oo. If we choose a step function for / so that m = we have J(l — cos 2 k 7rx)m(dx) = 

^°5Lq(1 — cos |j) < +oo so that the law of Nf is continuous and not Rajchman. 

4 Regularity results for multiple Poisson integrals. 

Let us first recall some links of our study with the Fock space. 
4.1 Random Poisson measure and Fock space. 

We recall that v is continuous (i.e. diffuse). Let us call simple the measurable functions / defined 
on {X m , X® m ) which are symmetric, finite sums of weighted indicator functions of sets of the 
form A\ x • • • x A m with disjoint A^s. 
On simple functions if we define 

Im{f)= / f{xi,...,x m )N(dxi)---N(dx m ) 

it is easily seen that 

^[Im{f)In{9)\ = Sm, n nl(j <?) L 2 (X™ ,X® m ,v*™)- 

Thanks to this equality I m (f) may be extended to / € L 2 (X m , X® m , u xm ) so that denoting / 
the symmetrized /, I m (f) = I m (f) and 

^[Im{f)In{g)] = 5 m ,nnl(f,g)L2(X™,X® m ,vX m )- 

Let us observe that for / G L 2 (X m , X® m , v* m ) the formula 

I m (f) = / f(xi, • • • , x m )l {Wif L j:Xi7 L Xj} N(dxi) ■ ■ ■ N(dx m ). 
1 There is an obvious misprint in this paper p!28 line 20 where = has to be change into 
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is a symbolic notation, because on the right hand side, the quantities to be substracted to the 
integral on X m are generally not defined for non regular functions /. 

It has a sense if / is well defined on diagonals by continuity, X being supposed topological. 
A sense may also be yielded by Hilbertian methods, supposing / allows to define trace operators. 

The sub-vector space of L 2 (VL, A, P) generated by the variables I n (f), f € L 2 (X n , X® n , u xn ) 
is the Poisson chaos of order n denoted C n . The equality 

L 2 (n,A,F)=R®t= 1 C n . (29) 

has been proved by K. Ito (see |21| ) in 1956. This proof is based on the fact that the set 
{N(Et) ■ ■ ■ N(E k ), (Ei) disjoint sets in X} is total in L 2 (tt,A,F). 

There are now several proofs of this result. A combinatorial proof is possible by counting the 
role of successive diagonals (cf |39| and )12| §4.1.) By transportation of structure, the density of 
the chaos has a short proof using stochastic calculus for the Poisson process on R + (cf Dellacherie- 
Maisonneuve-Meyer |17| p207). 

Thanks to the density of the chaos the following expansion is easily obtained (cf |45j ) for 
u £ L 1 n L°°{y) with small ||w||oo, 

+oo _. 

e iV(iog(i+«))-K«) = i + y -Uu® n ). (30) 

71=1 

Let us mention the relationship between the strongly continuous semigroup of the bottom struc- 
ture pt in L 2 (y) and the one of the upper structure Pi in L 2 (F) (see [T2] for a proof). For all u 
measurable function with — | ^ u ^ 0, 

Vt ^ 0, P t [ e N( - l °z( 1+u »] = e tf(iog(i-hptu))_ (3i) 

By (|30p and (|3 1 [) the vector spaces C n are preserved by Pt and 

Pt(I n (u® n )) = I n ((p t uf n )). (32) 

It is generally spoken of second quantization for the transform (pt) i— > (Pt)- More precisely the 
second quantization maps the generator a of pt to an operator on the Fock space which may 
be then lifted up either on the Wiener space or on the Poisson space and in this later case 
corresponds to the generator A of Pt. 

Remark 14. Let us suppose that the bottom semigroup pt be generated by a transition kernel 
Pt(x, dy) from (X, X) into itself, which be simulatable in the sense that there exists a probability 
space - that we choose here for the sake of simplicity of notation to be (R, 1Z, p) - and a family 
of random variables rjt(x,r) such that the law of ijt(x,r) under p(dr) be pt(x,dy) . 

Then, using our notation in which we have uj = J e x N(dx), the fact that the upper semigroup 
represents the evolution of independent particles each governed by pt and with initial law N (see 
the introduction of |12] ) may be expressed, for F measurable and bounded, by the formula 

P t F = EF( J £ m ( x>r ) N p{dxdr)) (33) 

in analogy with the Mehler formula for the Ornstein-Uhlenbeck semigroup on the Wiener space 
or extensions of it (see |7] pi 16). Applying f)33|) to F = exp AHog(l + g) for — \ ^ g ^ gives 



P t F = Eexp J log(l +g(rj t (x,r)) N Q p(dxdr) 
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what by formula (jl 3[) leads anew to (|3ip by a different way : 

P t F = expjVlog(/(l + g(rj t (x,r))p(dr)) = exp N log(l + p 4 c/). □ 

Remark 15. Surgailis |45] has shown that in the correspondence between pt and Pt given by 
( I32p a necessary and sufficient condition Pt be Markov is that pt and its adjoint be Markov 
operators (i.e. positivity preserving and s.t. ptl ^ 1). 

In our framework pt is selfadjoint and so is Pt. □ 

4.2 Decomposition of D in chaos. 

Let us precise some notation. On the upper space (fi, A, P, D, T) the Dirichlet form is denoted 
8. 

The product structure (X, X, v, d, 7 ) n will be denoted (X n , i/ xn , d n , 7n ) (cf [EE] Chap 
V). It is endowed with the Dirichlet form e n [f] = \ f -y n [f]dv. The functions in d n which are 
symmetric define a sub-structure of (X n , X® n , v xn , d n , j n ) denoted (X n , X® 1 ^, u xn , d n ^ sym , j n ). 
The semigroup associated with e n is denoted pf n . Our choice of gradient for the bottom space 
(see §2.3 above) induces a gradient for (X n , X® n , u xn , d n , 7„) that we denote (-)^ n with values 
in (L 2 (R,1l,p))® n : 

{f K ){x 1 ,r 1 ,x 2 ,r 2 , ■ ■ ■ ,x n ,r n ) = (f(-,x 2 ,-- ■ , x n ))\xi, n) + (f(x 1 , -,x 3 ,-- ■ , x n )) b (x 2 , r 2 ) H 

let us note that if / is symmetric, then f^ n is symmetric of the pairs (xj,rj). 

Let be f(xi,...,x m ) = /i(xi) • • • f m (x m ) G d m and 3(^1, . . . , x n ) = gi{x{) ■ ■ ■ g n (x n ) € d n . By 
polarization of ( j32j) Ptl m f = Impf m f gives 

£[im/,Jn0] = \{I m f ~ Ptlmf Jn9) L ^) = \{I m {f ~ pf™ f)Jn9) 

By the theory of symmetric strongly continuous contraction semigroups, we have F € B if and 
only if lim^o t £t[F] < +00 and £[F] = lim^m £t[^1- Taking / = g, we obtain that I m f £ D and 
£[im/] — ?Ti!e m [/]. Then by density we obtain 

Proposition 16. For f € d m the random variable I m f (= I m (f)) belongs to D. The vector 
spaces D m generated by I m f for f G d m , are closed and orthogonal in D. The sum 

ra> 1 

is direct in the sense of the Hilbert structure of IS) (|| • = || • ||? 2 + £[■])■ 
Every function F inD decomposes uniquely 

F = E[F}+ ]T I n (F n ) 

n > 1 

with F n € d n . 
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Proof. It remains only to prove the density of the Dirichlet chaos D n . Let be F £ D and let 
F = ^2 I n (F n ) be its L 2 -chaos expansion. Then 

\{F - p tF, F) L 2^ = \^Z n ^ 1 {In{F n -pf n F n ),I n F n ) L 2^ 

Since on the left-hand side t (F — PtF, F) ^ £[F] < +00 it follows that all terms on the right-hand 
side, which are increasing, possess limits what yields F n £ d n and the proposition follows. □ 

Let us emphasize that this proof is only based on the relation of second quantization (j32]) and 
would be still valid on the Wiener space for instance equipped with a generalized Mehler type 
structure (cf e.g. [7j pll3 et seq.) or on the Poisson space equipped with a non local Dirichlet 
form on the bottom space. 

Let u S L°° n d, applying the gradient operator (J to the two sides of (|30p gives 



JVlog(l+tu)-ti/(ti) 

n ^ 1 

what yields, taking terms in t n on both sides 

n-l 



(/„(«® n ))« = yV-l)« 7 rr/n-i-^"- 1 ^) / diV p. (34) 

(n-l-g)! 7 



and 



(35) 

If / is the symmetrized of fi(xi) ■ ■ ■ f m ( x m) then (|34p writes 

I m (ff = j (m/ m _ 1 / b -m(m-l)/ m „ 2 / b + m(m-l)(m-2)/ m „ 3 / b ----) dN p (36) 

where l m -p acts on the m — p first arguments of / and b acts on the last one, all free arguments 
being taken on the same point x. 

Extending formulae (|34[) - fj36 [ ) from tensor products to general functions / 6 d m supposes a 
priori that / does possess traces on diagonals. Indeed let us suppose / and g be regular so that 
values on diagonals make sense, then defining for regular symmetric functions f(x\, . . . , x m ) and 

7 

g(yi, ■ ■ ■ ,y n ) the (£;, ^^-contraction, for 1 ^ k ^ m and 1 ^ t ^ n, denoted / « g as follows 

7 

/ x 5 (xi, • • • ,x m _ fc ,yi, • • • ,y n -e,x) = 
k,£ 

j[f(xu--- ,x. m _ fc ,x,--- ,x,-),g{yi,- ■ ■ ,y n _ e ,x,--- ,x,-)]{x), 
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the function / x g is symmetric in {x\ } ■ ■ ■ ,x m _k) and in (yi, • • • ,y n _(). Then formulae (j34]) - 
(|36p extend to symmetric functions / and g as 



Wm(f), In®] = E T,(- 1 ) k+t ( rn _™Z-W Im - kIn - t / (/ X 9) (37) 
k=l£=l k ' 1 



where I m _& operates on the Xi's, I n -t operates on the yj's and N on x. 

But this formula is unsatisfactory because we know that I m (f) is defined and in D for general 
functions / E d m which dont have defined values on diagonals in general. Actually the values on 
diagonals cancel in formula (137j) . To see this we have to consider the Fock space for the gradient 
and to come back to the lent particle formula. 

The random Poisson measure NQp (cf §2.2) is defined on (£1 x f2, AqA, P X P) with intensity 
v x p on (X X R,X ® 1Z). It possesses an expansion in chaos : V-F £ L 2 (P x P) 



where J n denotes the multiple integral for NQp and where F n € L 2 sym ({v x p) xn ). 

Let us remark that the random Poisson measure N may be seen as a function of iV p and 
that the multiple integrals I n are nothing else but J n applied to a function G(x±,ri, ■ ■ ■ ,x n ,r n ) 
not depending on the r^'s. We can now state 

Proposition 17. Let be f 6 d m ^ ym , by Prop [16] i/ie multiple integral I m (f) belongs to D. 
Its gradient is given by 

Vm(f)f = J £-(/ m -i(/)) b ^ @p = mj I m -x(y) N p(dxdr) (38) 

where we note t/j(xi, . . . , x m _i, x, r) = (f(x±, . . . , x m -i, ■)) (x, r) and tf is defined as 

(p(xx, . . . ,x m _i,x,r) = ifj(xi, . . . ,x m _i,x,r)l{ a . i7 . :r vi=i,...,m-i} 

so that (p(-, ■ ■ ■ ,-,x,r) € L 2 J/m (z^ x ^ ?n ~ 1 )) and I m ^i((p) is defined, 
b) This gradient may also be written 

{Imff = JmiM (39) 

so that 

r[I m (f),I n (g)]=E[J m (f^)J n (g^)]. (40) 
Proof. Let be / € d.m,sym- Let us apply the lent particle formula to I m (f). We have 

e + I m (f) = I m (f) + mlm-tf P x z^-a.e. 
and since I-mif^ does not depend on x 

{e + I m (f)) b = m(I m _i/) b P x v x p-a.e. 
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Now, applying the operator e~ amounts to take the preceding relation with oj changed into e~u: 
and to work under the measure Ptv instead of P x v. That means that a functional F(N(u),x) 
is changed into 

e-(F(N(u),x)) =F(J u(y)l {y ^ x} N(dy),x) F N -a.e. 

Taking m = 2 for instance, we see that e~{I\f) must be written Pjv-a.e. J f(y,x)lr y ^t x \N(dy) 
instead of (I\f)(x) — f(x,x). Thus the part a) of the statement is a direct application of the 
lent particle formula. 

b) Since b takes its values in Lq(R, 1Z, p), it is equivalent to use the compensated random 
measure N Q p instead of N p in (j38j) . 

Now m J I m -i(ip)dN p = Jm(f^ m ) as seen by beginning with / = u® m , then polarizing to 
/ symmetrized of u% ® • • • ® u m and then to general / 6 d niS y m by density. □ 

Let us remark that formula ([38j) allows a new simple proof of the orthogonality of the chaos 
in D. Let / be as in the proposition. We have 

2S[I m f) = ET[I m f] = Em 2 J l[Im-l(fl{ Xi ^i})] N(dx) 
= m 2 j £- 7 [Im-l/] dNdF 
= m 2 J -f[I m _if]d¥diy (by Lemma [2]) 

= m 2 {m - l)\ j -i[f}dv x ^ m ~ 1 ^ dv 
= m\2e m [f]. 

and similarly with the scalar products. Now (|39p yields an even shorter proof using the orthog- 
onality of the chaos generated by J n under PxP, since (f m , g^ m ) L 2 (uxp) m = ^ e m[f,g}- 

Contrarily to the Wiener case the random variables I m (f) are n °t regular in general. Their 
distributions may contain Dirac masses. Even in the first chaos the fl or the T applied to 
I\u = Nu yields a non deterministic result, and the sharp operator does not diminish the order 
of the chaos. Studying regularity of multiple integrals needs therefore additional hypotheses. 

4.3 Density for {h(g), I n {g® n )). 
Relation fj34[) yields immediately 

Let us denote I the column vector of (Ji, . . . , I n ), we have 

r[z,x*] = J vvHa] dN 

with V the column vector of (-1, -I x + g, . . . , n! Y2=i J^iy.i- 1 ) 1 *^' 1 )- 

Let us precise now some hypotheses. We suppose z^{7[<7] > 0} = +oo and that assumptions are 
fulfilled such that we have (BC) on the bottom space and (EID) on the upper space, as usual. 
If for some uj € f2 the matrix F[I,I t ] is singular, this means that all the vectors 

V(u,Xi(ui)) for Xi € supp(w) n {j\g} > 0} 
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belong to the same hyperplan of M. n , in other words, this implies that there exist Xq(oj), . . . , A re _i(w) 
not all null such that: 

n T 

on all the points of supp(w) D {"fig] > 0}. 

Since g £ d, by (EID) on the bottom space — which is always true for scalar functions — the 
measure 5 f *[l{ 7 [ 9 ]>o}- zy ] is absolutely continuous hence continuous (diffuse). As f{7[<?] > 0} = +00 
the random Poisson measure image by g of the points of N which are in {"f[g] > 0} do possess 
infinitely many distinct points. Hence the g(Xi(iv)) cannot annul a polynomial except if it is 
identically sero. 

The question is therefore to know whether 

n 

-AoM + Ai(w)(-/i + x) + ■ ■ ■ + A n _i(a;)n! V ^-i = 

^— ' [n — K)l 

implies Xo(ui) = • • • = A n _i(u;) = 0. 

But this is due to the fact that the annulation of the coefficients of this polynomial builds a 
triangular linear system whose diagonal terms are — Ao(w), . . . , n!( — l) n A n _i(w). We have proved 

Proposition 18. If the upper structure satisfies (EID), for g £ L°° R d such that 
u il[g] > 0} = +00 the vector (h(g), ■ ■ ■ ,I n (g® n )) has a density on W 1 . 

Remark 19. This result is quite different from what happens on the Wiener space since there 
the law of (I\(f), . . . , I n {f® n )) is carried by the algebraic curve of equation 

' x 2 = 2\H 2 {\\f\\\x 1 ) 

k x n = n!i? n (|[/|[ 2 ,3;i) 
where H n (X,x) is the Hermite polynomial given by 

t 2 X 00 
exp(te -—) = J2 t n H n (X, x). 

n=0 

4.4 Density for (I ni {f?*), . . . , I^ff*)). 

Let /=(/!,..., / p ) € (L 1 n L°° n df and let be J the column vector (J m (/f ni ),..., I Up (f® np ) ) 
where we suppose 1 Vi. 

Defining the polynomials Vi by Vi(x) = i\ fX)fc=i [i-k)\ (~ l)*^ -1 ) we have by ([ITjl the 
equality between p x p-matrices 

T[J,J t } = j {Vn i (fi)'Pn j (fjhlfiJj}) ij dN. 

By Lemma [8] 

{detriJ,^] = 0} C { / det7[/,/'Fn 1 (/i) 2 ---^n p (/ P ) 2 dN = 0}. 
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Let us assume v{det^/[f,f t ] > 0} = +00, and that we have (EID) below and above. The image 
by / of l{det[/,/*]>0} ' v 1S absolutely continuous w.r. to Lebesgue measure and the Poisson random 
measure image of N |{dct [/,/*] >o} has an absolutely continuous and infinite intensity measure, it 
possesses necessarily points outside the finite union (less than ^|=i( n « — 1)) of hyperplans defined 
by V ni {xi) = whose term of highest degree is (rij)!(— \) ni x™ % ~ . We obtain 

Proposition 20. If (EID) holds below and above, and if j[f, /*] is invertible u-a.e. 
(I ni (ff' ni ), . . . , In p (fp np )) has a density as soon as ri{ ^ 1 Vi. 

Remark 21. Let us compare with the situation on the Wiener space. We dispose only of 
sufficient conditions of regularity, but we can nevertheless compare the thread of the arguments. 
We have DI n {g® n ) = n/„_i(^® (n-1 ))s and 

r^r*),^-^)] =n i n^„ i _ 1 (/f^- 1 ))/ nj ._ 1 (/ n ^f 1 )) f f ifj dt. 

Since (EID) holds on the Wiener space a sufficient condition of density of J is that almost surely 
the vector 

(niln.-lC/f^-'ViW, • • • > V nr i(f rl) )/pW) 

generates a p-dimensional space when t varies. It is easily seen by induction on n that 
V/ G L 2 {dt), 11/11 / : F{I n {f® n ) = 0} = 0. It follows that on the Wiener space, J has a 
density as soon as ni ^ 1 Vi and . . . , f p ) are linearily independent in L 2 (dt). 

4.5 Other functionals of Poisson integrals. 
Density of (JV(/ifo), . . . , N(f n (g))). 

Let 5 £ L°°nd and let fi be regular real functions on R. Let us denote /C = (N(f\ (g), . . . , N(f n (g))) t 
and suppose i/{ 7 [ ff ] > 0} = +00. From I^/ifo)), JV(/ite))] = f f^ffaWgjdN we obtain 
that the matrix T[/C, /C*] is singular if the vectors (/((<?), ■ ■ ■ , f ' n (g)) taken on the points of u are 
in a same hyperplan. Now the points g(x), x 6 supp(w), have an accumulation point at zero. 
We obtain 

Proposition 22. Suppose (EID) holds above, g € L°°f]d, ^{7 [5] > 0} = +oo, and the functions 
fi be analytic at the neighborhood of O such that (1, f%, . . . , f n ) be linearily independent, then 
(N(fi(g)), N(f n (g))) has a density. 

Since there are infinitely many distinct points g(x), x € supp(w), we see also that without 
analyticity hypothesis it suffices that any hyperplan cuts the curve (f[(t), . . . , /^(£))teK a t a 
finite number of points, the /j being supposed C 1 H Lip. 

Density of (£. N(fj), Ej( N (fj)) n )- 

Let us consider <I> the column vector of the polynomials 
<&fc(xi, . . . ,x n ) = X]j=i X ji' / = (/l)---i/n) ^ d" and let us pose V the column vector of 
the $ fc (iV(/i), . . . , N(f n )). We obtain 

Y[VX] = V1>(Nf 1 ,...,Nf n )T[Nf,Nf t ](V^) t (Nf 1 ,...,Nf n ) 
= V*(Nf u Nf n ) J 7 [/, /*]diV(V$)*(J\T/i, . . . , iV/ n ) 
detr[V,V*] = (detV$(iY/ 1 ,...,iY/„)) 2 det/ 7 [/,/^A r 
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where V$ is the Jacobian matrix of <3?. 

det V3> is a Vandermonde determinant, if v{fi ^ fj} = +00, det V$(A r /i, . . . ,Nf n ) cannot 
vanish by Paul Levy's theorem. 

J 7[/j f l \dN is an infinite sum of non negative symmetric matrices, as before we can state 

Proposition 23. Supposing (EID) above, u{fi 7^ fj} = +00 \fi 7^ j, and z^{det 7 [/, /*] > 0} = 
+00, then V has a density. 

4.6 Density of /„(/) for / e d 

If f(x±, ■ ■ ■ ,x n ) is a symmetric element of d n , the function (f(x\,--- , x n _i, •)) (x,r) may be 
seen as a symmetric Hilbert valued function in d n _i(i?) with H = L 2 (u x p). So that we can 
iterate the operator b going down on the arguments (f(x\, • • • , x n _2, ■, ■))'''' S d n _2{H <S> H). 
Let us apply this with Prop [TF] 

T[I n (f )] = n 2 J e- 7 [J„-i(/)]dtf. 
By Lemma [2]for r[/ n (/)] to be > it suffices 

7 [/ n -i(/)] > Fx z/-a.e. 

i.e. that l n -\f^ be 7^ P x v x p-a.e. hence it suffices that f x /5-a.e. I n -iy have a continuous 
law. 

Now getting down the induction and using Paul Levy's theorem yields that it suffices that 

{y x p) n - x -a.e. v{xi : / {n " 1)b (xi, x 2 , r 2 , . . . , x n , r n ) / 0} = +00. 

Applying this to the classical case where the bottom space is equipped with the Lebesgue 
measure and the form e[f] = | J f' 2 (t)dt, where we can choose / = /' • £ with £ reduced 
Gaussian, we obtain 

Proposition 24. For n ^ 2, I n {f) has a density if the Lebesgue measure of the set {x\ : 
Qx 2 ...J x 7^ 0} i- s infinite dx-i ■ ■ ■ dx n -a.e. 

This extends to the classical case on M. d taking f^ = + • • • + ^-£,n with the £j i.i.d. 

reduced Gaussian. 

Remark 25. There is a major difference with the case of the Brownian motion about the sum 
of the series 

In the case of Wiener space this sum is a function of f fdB = I\{f) since it is equal to 
e tJfdB-±t 2 \\f\\\ Qn the Poisson space it is not a function of h(f) = N(f) but of iV(log(l + tf)) 
and for / G L°° n d and small t by our usual argument using Paul Levy's theorem the pair 
(Nf,N log(l + i/)) do have a density if i/{7[/] > 0} = +00. 

It is natural to ask about the density of the vector (iVlog(l + t\f), . . . , AT log (1 + t n f)). For 
/ G L°° n d, supposing < t\, . . . , t n < ||/||ooi by the method it suffices to have (BC) down, 
(EID) above, ^{ 7 [/] > 0} = +00 and the t{ to be distinct. 
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Remark 26. In the Wiener case multiple integrals obey a product formula (cf. Shigekawa |43] 
p276) allowing to express explicitely / m [/]/ n [g(] as linear combination of multiple integrals of 
order less or equal to m + n. 

A similar formula exists on the Poisson space slightly more complicated. It may be obtained 
in the following way. Let u,v € L 2 n L°°(v) with small uniform norm. By the relation 

JV(log(l+su))— sv(u) iV(log(l+to)) — tv(v) (\og(l+ su+tv+stuv)) — v((su+tv+stuv) stv(uv) 

thanks to f|30[) we have 

(1 + V ^ T / m (u® m )(l + V -I n (v® n ) = (1 + T -J P ((su + tv + stuv)®ne stu{uv) 

m=l n=l p=l 

and the product formula is obtained by identification of the term in s m t n of the two sides. Then 
it may be extended by polarization to / and g f or / = f± ® • • • ® f m and g = g\® ■ • • ® g n an d 
then for general / € L 2 (z^ xm ), g € L 2 (v xn ) by density. See [22] , [35] , |46] for different forms of 
such a formula, also [40], [39], and [T7] p261 for a general expression and proof. 

If we apply this product formula to J m (f' 9rn )Jn(g^ n ) using EJfc(fo) = hp(dr\) ■ ■ ■ p{drk)) 
for h(xx,ri, ■ ■ ■ ,x^,r^) € l? sym (y x /9) xfc ) we could obtain another expression of T[I rn f, I n g] = 

EJ m (f bm )J n (g K ) to be compared with fl40j) . 
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